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ABSTRACT 

^  Let  L,  R,  Y  be  arbitrary  real  constants.  A  wing  with  fixed  span  which 
produces  lift  L,  roll  moment  R,  yaw  moment  Y  and  which  has  minimal  induced 
drag  D  is  wanted.  This  problem  arises  in  airplane  engineering.  It  is  solved 
by  means  of  Prandtl's  lifting  line  theory  combined  with  computation. 
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SIGNIFICANCE  AND  EXPIANATION 


"Adverse  yaw"  has  been  with  us  since  the  beginning  of  the  era  of  heavier- 
than-air  flight.  The  pioneers  Lilienthal ,  the  Wright  brothers  and  others  must 
have  encountered  this  phenomenon.  It  consists  of  the  fact  that  an  airplane 
Initiating  a  turn  normally  moves  in  the  wrong  (adverse)  direction  if  only 
stick  movement  and  no  rudder  deflection  is  used. 

It  appears  that  this  phenomenon  has  not  been  studied  very  much  by 
aerodynamicists  or  airplane  engineers.  There  are  several  reasons  behind  this 
fact:  The  first  may  be  that  adverse  yaw  is  normally  more  of  a  nuisance  than  a 
real  problem.  In  most  airplanes  it  can  be  controlled  with  the  rudder.  When 
teaching  people  how  to  fly,  it  is  standard  to  demonstrate  the  adverse  yaw  and 
its  manual  correction. 

Another  reason  may  be  that  adverse  yaw  is  most  powerful  for  wings  with 
high  aspect  ratio,  e.g.  sailplanes  and  airliners.  Sailplane  manufacturers, 
however,  usually  don't  invest  in  scientific  research.  While  the  airline 
industry  could  make  such  an  investment,  it  has  not  done  so  because  the  auto¬ 
pilot  takes  care  of  adverse  yaw  during  flight. 

Nevertheless,  in  both  situations  this  correction  for  adverse  yaw  expends 
unneccessary  energy.  It  is  the  aim  of  this  paper  to  show  how  to  reduce  the 
energy  required  in  the  process  of  correcting  for  adverse  yaw  by  minimizing  the 
resulting  additional  drag. 

The  interest  of  the  author  in  this  subject  comes  from  his  work  with 
tailless  airplanes  where  adverse  yaw  is  a  serious  problem.  It  is  hoped  that 
other  kinds  of  airplanes  may  also  benefit  from  the  results  presented  in  this 
report . 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


MINIMAL  DRAG  FOR  MINGS  WITH  PRESCRIBED  LIFT,  ROLL  MOMENT  AND  YAH  MOMENT 


Karl  L.E.  Nickel 


1.  Explanation  of  the  physical  problem 


Every  airplane  pilot  ia  faailiar  with  the  following  phenomenon:  Suppose 
the  pilot  wants  the  airplane  to  roll.  Then  by  using  the  ailerons  as  sketched 
in  Figure  1  a  roll  moewnt  R  is  produced.  However ,  more  than  the  (desirad) 
roll  around  the  longitudinal  axis  occurs.  The  airplane  also  turns  its  nose  to 
the  side  as  if  a  rudder  deflection  had  been  made,  see  Figure  2.  This  is  known 
as  "adverse  yaw".  In  most  cases  it  is  not  wanted  and  is  regarded  as  a 
nuisance. 
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Figure  1 .  Sketch  of  airplane  with  lift  L  and  roll  moment  R  produced  by  aileron 
deflection. 

! 


Figure  2.  Sketch  of  airplane  with  yaw  moment  Y  around  the  vertical  axis 
produced  either  by  the  rudder  deflection  shown  or  by  the 
aileron  deflection  shown. 

Sponsored  by  the  United  States  Army  under  Contract  No.  FaaG29-80-C-0041 . 


The  physical  reason  for  adverse  yaw  la  that  on  the  side  with  the  aileron 
bent  downward  not  only  the  (local)  lift  is  Increased,  hot  that  there  la  also 
an  additional  (local)  drag.  This  drag  is  called  (by  Prandtl)  the  'induced* 
drag.  It  pulls  this  part  of  the  wing  back  and  therefore  produces  the  yaw 
moment  Y. 

From  this  discussion,  the  following  property  is  plausible:  The  adverse 
yaw  is  proportional  to  the  lift  L  produced  by  the  wing.  Hence,  no  adverse  yaw 
occurs  if  L  is  sero  (vertical  dive)  and  if  the  aileron  configuration  is 
symmetric.  For  negative  lift  L  (upside-down  f lying)  strong  downward  gust) 
but  unchanged  roll  moment  R,  the  sign  of  the  yaw  moment  Y  is  obviously 
reversed. 

In  most  airplanes  of  today  (or  of  the  past)  the  adverse  yaw  is 
annihilated  by  the  (auto-)  pilot  just  by  steering  against  it  with  the  rudder. 
There  exist,  however,  at  least  two  kinds  of  airplanes  at  present  for  which 
this  is  difficult  to  do  or  even  impossible: 

1)  Tailless  airplanes  with  no  vertical  tail  fins  and  with  no  rudder. 

Most  of  today's  hanggliders  belong  to  that  category. 

ii)  High  efflcency  sailplanes  with  their  usually  very  large  wing  span 
(24.5  m  or  80  ft  in  1983).  In  some  of  these  airplanes,  at  maximum 
lift  full  rudder  deflection  is  not  sufficient  ( I )  to  fight  adverse 
yaw.  Hence,  even  though  adverse  yaw  is  a  very  old  problem,  it  seems 
to  have  some  modern  significance. 

There  are  several  constructive  devices  to  diminish  or  to  annihilate 


adverse  yaw.  They  will  not  be  described  here.  Because  of  the  properties 


mentioned  above  they  are,  however,  normally  effective  only  in  a  certain 
limited  range  of  the  lift  L.  - 

Since  adverse  yaw  ia  generated  by  an  asymatric  distribution  of  the  local 
induced  drag  one  may  suspect  that  it  could,  perhaps,  be  nullified  by  an 
adequate  modification  of  that  distribution.  This  leads  to  the 

First  Question  :  Are  there  special  kinds  of  ailerons  and/or  wings  for 
which  adverse  yaw  does  not  exist  for  any  values  of  lift  L  and  roll 
mosient  R  ? 

It  will  be  shown  in  this  paper  that  there  are  infinitely  many  such 
configurations.  This  is  a  very  'clean*  solution  to  the  problem  of  fighting 
adverse  yaw  since  no  additional  devices  are  needed:  The  wing  itself  and  its 
ailerons  will  be  shown  to  suffice  for  the  elimination  of  adverse  yaw. 

Obviously  this  solution  will  not  be  'for  free'.  In  order  to 
counterbalance  the  one-sided  induced  drag  one  has  to  add  some  additional  drag 
on  the  other  side  of  the  wing.  Hence  the  total  induced  drag  D  of  the  airplane 
is  increased.  This  leads  to  the 

Next  Question  :  Is  it  possible  to  select  one  solution  to  the  First 
Question  which  produces  minimal  total  drag  D  ? 

It  will  be  shown  that  there  is  exactly  one  such  solution. 

In  most  cases  an  airplane  rolls  if  it  wants  to  take  a  turn.  In  this  case 
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roll  and  positive  yaw  are  wanted  at  the  sane  time.  Hence  one  may  ask  if  one 
can  replace  rudder  deflection  <  full  or  partial  )  by  a  favorable  distribution 
of  the  induced  drag  on  the  wing.  In  extending  and  generalizing  these 
Questions  one  cones  therefore  to  the 

Final  Question  i  Are  there  configurations  of  wing  and  ailerons  with 

minimal  drag  D  for  arbitrarily  prescribed  values  of  lift  L  ,  roll  moment 

R  and  yaw  moment  Y  ? 

It  will  be  shown  that  this  question  has  in  general  (  i.e.  for  R  +  0  )  exactly 
one  solution.  Examples  will  be  presented.  Moreover,  explicit  formulae  will 
he  given  which  bound  the  minimal  drag  from  both  sides  within  an  error  of  less 
than  2.1%  . 

The  approach  used  in  this  paper  will  be  Prandtl’s  lifting  line  theory 
(see  [4] ) •  This  theory  is  a  first  approximation  to  the  lifting  surface  theory 
and  normally  cannot  be  used  for  general  wing  shapes,  e.g.  swept  wings.  For 
the  prohlem  considered  here,  however,  it  can  be  used  without  restriction  due 
to  Hunks's  shifting  theorem  [1].  This  states  that  for  the  evaluation  of  the 
total  Induced  drag  only  the  projection  of  the  lift  distribution  over  the  wing 
onto  the  Trefftz-plane  is  essential,  regardless  of  sweep. 

For  the  same  reason,  only  the  lift  or  circulation  distribution  over  the 
wing  is  considered  in  what  follows.  There  is  (  for  R  f  0  )  exactly  one  such 
optimal  distribution.  There  are,  however,  infinitely  many  wing  and  aileron 
combinations  which  produce  this  desired  distribution.  The  connection  between 
wing  and  aileron  shapes  and  circulation  distribution  is  known  and  can  be 
handled  by  standard  methods. 

Also  not  treated  here  is  the  following  question  of  optimality:  By  using 
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the  final  results  of  thla  paper ,  one  can  design  an  airplane  with  a  vary  small 
rudder.  Thla  reduces  drag  In  straight-  ahead  flight.  During  turns,  however, 
a  larger  drag  is  produced  due  to  the  generation  of  positive  yaw.  What  is  the 
best  combination  of  favorite  yaw  and  small  rudder  ? 

The  results  of  the  forthcoming  theory  have  been  field-tested.  This  was 
done  by  constructing,  building  and  flying  a  specially  designed  ultralight 
airplane  called  ’Palter  1*  (Palter  means  ’butterfly*  in  German  but  is  also 
used  for  ’foldable*  -  the  pun  was  intended).  The  Palter  1  is  a  tailless 
aircraft  with  swept  wings  which  has  no  vertical  fins.  It  uses  only  ailerons 
for  the  steering  around  the  three  axes.  It  is  quite  obvious  that  both  roll 
moment  and  pitch  moment  can  be  generated  by  using  ailerons  (  the  latter 
because  of  the  swept  back  wings  ).  But,  on  the  other  hand,  it  looke  at  first 
as  if  the  yaw  moment  (which  is  necessary  for  turns)  cannot  be  produced  in  this 
way.  In  the  past,  all  known  flying  wings  which  were  steered  by  aerodynamic 
means  (that  is,  not  by  shifting  weight,  as  the  hanggliders  do  it)  used  either 
vertical  fins  or  drag  rudders  for  producing  the  yaw  moment  (  Lippi sch,  Horten, 
Northrop  et  al.)i  not  so  with  the  Palter  1.  It  is  steered  around  the 
vertical  axis  with  only  the  ailerons  which  produce  the  desired  yaw  moment.  By 
using  the  theory  derived  below,  care  was  also  taken  to  make  the  (unavoidable) 
additional  induced  drag  as  small  as  possible. 
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2.  Th»  Math— tlcal  grobli 


2.1  Definition* 

The  following  notations  and  definitions  are  standard: 

b  ”  span  of  the  wing, 

S  *  surface  of  the  wing, 

2 

X  :■  b  /S  »  aspect  ratio, 
v  ”  velocity  of  the  aircraft, 
p  ■  air  density, 

q  :■  »  kinetic  energy,  soaetimes  called  'dynamic  pressure*, 

L  -  lift, 

D  “  (induced)  drag, 

R  ”  roll  moment, 

Y  »  (induced)  yaw  moment, 

CL  i»  L/qS, 

CD  «“  D/qS, 
b 

CR  ^*2  ' 
cy  s»  Y/qs|  , 

x,y  ■  coordinates  in  wing  direction.  In  what  follows  these  coordinates 
are  (by  an  affine  transformation)  restricted  to  the  interval 
-1  <  x,y  £  +1 . 
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2.2  Prandtl's  Lifting  Line  Theor1 


Definitions  And  Properties 


Z(y)  -  (local)  circulation  density 

-  (local)  lift  density  at  the  wing. 


at<y) 


+1 

1  X  dx 

2  dx  y  -  x 


1 


) 


-  (local)  induced  angle  of  attack. 


With  these  two  definitions  one  gets  for  the  above  defined  coefficients  the 
relations 

+  1 

cL  -  X  J  Z(y)  dy  , 

-1 

+  1 

c„  -  X  /  Z(y)  ydy  ,  2) 

-1 

♦  1 

cD  -  X  /  Z(y)  (^(y)  dy  , 

+  1 

Cy  -  X  /  Z(y)  a^y)  y  dy  .  2) 


2 . 3  Transformation  To  An  algebraic  Problem 


With  the  transformation 

y  ”  cos  s  for  0  <  s  <  s 


’)  As  usual,  the  letter  c  Inside  an  integral  sign  means  the  cauchy 
principal  value  of  that  Integral. 

2)  Depending  upon  the  orientation  of  the  axis  system  used  sometimes  in  the 
literature  cR  and  Cy  are  replaced  by  -cR  and  by  -Cy  . 


-  7  - 


and  by  developing  Z  Into  a  Fourier  aariaa 


Z(coa  a)  •  2  £  a  Bin  va 


on*  gets 


cL  -  »Aa1  . 


=R  -  IX*2  * 


V»1 


*x  J 


v-1 


■v  ' 


TX  I  <2v+1>«v*vt,  • 


V«1 


In  order  to  get  rid  of  the  constants  rX  and  — X  ,  it  is  useful  to  define  new 
dimensionless  forces  and  moments  by 


L’  L/wXqS  -  c  /*X  , 

L 

R’  l-  R/^XqSj  -  2cj/’X  ' 
O'  D/wXqS  -  CpAX  * 
Y’  i-  Y/^Xqs|  -  Zc^nX  . 


-  8  - 


[ 


Than 


R'  “  «2 • 

m 

D*  *  I  u«v2  * 

V-1 

oo 

Y*  -  l  (2v+1)ayav+  . 

v-1 

In  what  follows  the  abbreviation  a  -  {a^}  for  the  sequence  {a^}  will 
sometimes  be  used;  moreover,  the  reduced  drag  D*  will  be  regarded  as  a 
functional  D'(a)  acting  on  a,  ms  leads  to  the  definition  of  the  set  A  of 
admissible  sequences  by 

A  i-  {a | D' (a)  <  “}  . 

The  Cauchy-Schwarz  Inequality  Implies 

OD 

I  £  (2v+1)avau+1 1  <  <®  for  any  a  e  A  . 

v-1 

with  these  definitions,  the  problem  posed  in  the  Final  Question  of 
Section  1  becomes  : 
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Probigg  i  Given  are  the  real  number a  L' ,  R* ,  ¥’. 

Wanted  la  an  admissible  sequence  a  £  A  which  aatlafiea  the 
conditiona 

(1)  a1  -  L' , 

(2)  a2  -  R' , 

«P 

(3)  l  (2v+1 )a  a  -  Y* 

*-  V  v+1 

v-1 

such  that 

oo 

p  2 

(4)  D'(a)  2.  va  «  minimum  . 

w-1 

The  existence,  and  (for  R'  f  0)  the  uniqueness,  of  a  solution  a  e  A  to  this 
problem  will  be  proven  in  Section  5. 
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3.  Some  Preliminary  Results 

The  results  presented  In  this  section  are  all  quite  simple.  Surprisingly 
enough,  the  author  did  not  find  them  in  papers  or  books  on  aerodynamics  or 
airplane  engineering  (with  the  exception  of  Sections  3.1  and  3.2).  It  may, 
therefore,  be  worthwlle  to  print  them  here. 

3 . 1  Only  The  Lift  L  Is  Prescribed. 

In  this  case  the  solution  a  to  the  above  problem  is  obviously 

al  ’  L'<  av  *  0  for  u  “  2<  3 . 

This  problem  and  its  solution  goes  back  to  Prandtl  [4]  in  1918.  In  physical 
terms  it  reads  as 

Theorem  1  (L. Prandtl):  The  optimal  circulation  distribution  on  a  wing 
with  prescribed  lift  L  is  elliptic  and  is  given  by 

2cl  / - T 

(5)  Z(y)  77  /  1  -  y*  • 

The  coefficient  of  the  induced  drag  for  this  case  is 


2 


3.2  Both  Lift  L  And  Roll  Moment  R  Are  Prescribed,  But  Nothing  Else. 

This  is  a  special  case  of  a  more  general  problem,  solved  by  the  author  in 
1951,  see  [3).  The  solution  a  to  this  problem  can  also  be  found  quite 
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easily.  It  Is 


»1  "  L**  “2  *  R'<  «y  -  0  f°r  V  “  3'  *' 


Expressed  in  physical  terns  it  is  stated  as 


Theorem  2  (K. Nickel):  The  optimal  circulation  distribution  on  a  wing  with 
both  prescribed  lift  L  and  roll  moment  R  is 


2c. 


(7) 


Z1.2(*’ 


-jf  /7 

TC  X 


— -  8c 

2  .  _ R 

»  +  77 


/TT7 


It  is  a  linear  combination  of  a  (half)  ellipse  and  a  (half)  lemniscate. 
The  coefficient  of  the  induced  drag  for  this  optimal  distribution  Is 


(8) 


i  X  * 


Obviously  the  problem  solved  in  Section  3.1  is  a  special  case  of  this  result 
if  cR  “  0.  The  adverse  yaw  which  is  produced  by  this  distribution  can  be 
found  from  equation  (3).  Its  coefficient  cy  is 

CY  "  77  CL  CR  * 

It  shows  that  -  at  least  for  this  case  -  the  adverse  yaw  is  proportional  to 
both  lift  L  and  roll  moment  R,  as  remarked  above. 
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3.3  Non-Optlmal  Distributions  for  Given  Lift  L  And  Roll  Moment  R  With 
Vanishing  Mvwm  Yaw  (Y  -  0). 

There  are,  obvioualy,  infinitely  many  sequences  a  c  A  for  which  Y  «  0, 
i.e.  with 


(9)  l  (2vt1)avav  -  0  . 

v— 1 

Only  two  of  thee  will  be  considered  here. 


3.3.1  ay  «  0  for  v  -  4,  5,  ...  . 


In  this  case  the  equation  (9)  becomes 
3  *!  *J  H  *2  *3  ’  0  • 

Hence,  by  assuming  a2  f  0,  w.l.o.g.  one  gets 
3 

*3  S"  "  5  “l  ' 

This  solves  (9)  for  any  value  of  a2>  1*«»  o f  the  roll  moeient  R.  The 
corresponding  circulation  distribution  is 


(1°)  (y) 


8cl  „  ,  2,  rr~i .  8cr  r~~2 

5wX  (2“3y  ’  7  1_y  +  7T  y  7  1_y 


The  induced  drag  coefficient  is  given  by 


(11) 


2  2 
°L  CR 

cd  ■  2-08  nr  ♦ 8 


w  X  ' 


For  R  «  0  it  exceeds  the  smallest  possible  induced  drag  for  a  given  lift 
(see(6))  by  the  factor  of  2.08  (I).  This  is  the  price  paid  for  annihilating 
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the  adverse  yaw.  It  will,  however,  be  shown  later  that  there  are  “better" 
distributions  for  which  there  is  no  additional  drag  for  R  •  0  .  A  very  simple 
example  of  this  is  shown  in  the  next  Section 

3.3.2  Sj  »  0  and  ay  »  0  for  v  «  6,7,  ...  . 

Free  the  equation  (9)  one  gets  the  condition 

*4  *5  ’  "  I  *1  a2  ' 

with  eoaie  freedom  in  choosing  a^  or  as  . 

By  putting 

•4  *  -“«5 

and  by  determining  a  so  that  the  induced  drag  becomes  minimal  one  gets 


Obviously  this  distribution  is  much  "better"  than  that  given  in  3.3.1  .  Tor 
R  -  0  (  no  roll  moment)  Its  drag  coefficient  reduces  to  that  value  given  in 
(8).  For  R  *  0  it  exceeds  the  optimal  value  in  (8)  only  by  the  term 
8  /5  |  cR  I  /  3  *  X  which  is  small  for  small  cL  and/or  cR  . 

The  corresponding  circulation  distribution  is 
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from  equation  (7) 


(13) 

It  la  sketched  in  Figure  3.  By  using  the  distribution 
in  section  3.2  together  with  Z45  froa  (13)  one  gets  the  very  inter  eating 

Result i  Both  circulation  distributions 
Z(y)  i-  zi.2(y>  ±  CL  cr  z4.5(y) 

have  the  same  lift  L,  roll  moment  R,  vanishing  yaw  moment  Y  ■  0  and 
induced  drag  as  given  in  (12). 

See  the  Figures  3  and  4. 


Figure  3.  The  circulation  distribution  Z^j. 


Figure  4.  the  three  circulation  distributions  z1>2  and  Z,>2  i.  z4.5  *  They  have 
all  the  same  lift  L  and  roll  moment  R.  The  latter  two  have  both 
vanishing  yaw  moment  V  *  0  and  attain  the  same  induced  drag  D. 
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In  Figure*  5  to  8  some  circulation  distribution*  are  presented  which 
solve  the  problesi  of  this  paper.  In  all  cases  Y  “  0,  i.e.  null  adverse  yaw  was 
chosen.  Because  the  problem  is  then  homogenuoua,  without  loss  of  generality 
L'  : ”  1  was  assured.  In  these  Figures  four  typical  roll  moorents  were 
selected.  These  results  have  been  computed  by  Mrs.  Norbert,  and  the  Figures 
were  drawn  by  Mrs.  Sturm,  both  Freiburg  i . Br . /GERMANY .  The  method  used  is 


In  Figure  S  the  four  optimal  circulation  distributions  are  shown.  In 
order  to  see  more  clearly  the  details,  these  distributions  are  broken  up  into 
the  corresponding  symmetric  part  Z0  and  the  antisymmetric  part  Z by  the 
equation 

Z  -  Zg  ♦  za  . 

These  two  parts  are  presented  in  Figures  6  and  7.  Ks  cm  bee  seen  from 
Figure  6,  the  optimal  antisymmetric  circulation  distribution  for  the  largest 
roll  moment  R'  *  2  looks  very  similar  to  the  lemniscate  which  is  the 
optimal  distribution  if  the  condition  Y  «  0  were  removed  (  see  Section  3.2). 
It  can  be  shown  that  this  is  typical,  i.e.  that 

I  2 

Z^ly)  — ■*  const,  y  /  1  -  y  as  R— ♦«. 


The  proof  of  this  will  be  given  elsewhere. 


Figure  6.  The  antisymmetric  part  of  the  distributions  of  Figure  5. 
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Flcrure  7.  The  symmetric  pert  of  the  distributions  of  Figure  5. 


In  order  to  get  a  closer  look  at  the  symmetric  circulation  distribution 
in  Figure  7  a  new  symmetric  distribution  with  sero  lift  is  defined  by 

2cl.  /  2 

2  „<y)  i-  2  (y)  -  ~~  /  1  -  y  . 

SO  S  »  A 

It  is  presented  in  Figure  8.  Mote  that  this  distribution  gives  negative  ( I ) 
additional  load  at  the  wlngtips> 


19 


Similarly  to  the  earlier  result  it  can  be  ahown  that  asymptotically 


Z(y)  -♦  Z1.2.3(y)  as  R  -*  » 

for  all  optimal  circulation  distributions!  where  is  defined  in 

equation  (10)  . 

3.5  Some  Simple  General  Results. 

Let  again 

z(y)  -  za(y)  +  za(y) 

with  Zg  and  Zg  the  symmetric  and  the  antisymmetric  parts  of  Z, 
respectively.  It  is  known  since  Prandtl  [4]  (see  Section  3.1)  that  the 
elliptical  symmetrical  distribution 

(5)  V*’  =■ 

is  optimal  with  respect  to  the  drag,  if  no  conditions  on  roll  moment  R  and  yaw 
moment  Y  are  required.  Hence,  one  may  ask  if,  given  the  symmetric  distribution 
in  (5),  an  additional  antisymmetric  distribution  Za  can  be  found  such  that 
Y  •  0,  l.e.  such  that  no  adverse  yaw  occurs.  By  looking  at  (9)  and  bearing  in 
mind  that  the  Fourier  coefficients  aJv  belong  to  the  antisymmetric  and  the 
*2v+1  t0  symmetric  distribution  one  sees  immediately  that  no  such 
distribution  Zg  can  be  found,  provided  that  L  +  0  and  R  +  0.  Hence,  one  has 
the 

Theorem  3s  Assume  L  +  0,  R  +  0.  There  is  £0  circulation  distribution  Z(y) 
with  the  symmetrical  part  Zg  as  In  equation  (5)  for  which  no  adverse 
yaw  occurs. 
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Now  a  somehow  'dual*  question  arises  i  It  is  known  (see  Section  3.2) 
that  the  leaniscats  shaped  antisymmetric  distribution 

8  CR  / - 2 

(14)  za(y)  i-  — y  /  i-y 

is  optimal  if  there  is  no  condition  on  the  yaw  moaMnt.  In  analogy  to  the 
above  one  aay  ask,  therefore,  if  there  are  symmetric  distributions  Z( 
corresponding  to  Za  frost  (14)  such  that  Y  *  0  .  By  looking  at  equation  (9) 
one  sees  easily  -  opposite  to  the  above  result  in  Theorem  3  -  that  there  are 
infinitely  many  such  distributions.  Che  can,  therefore,  try  to  solve  the 
basic  problem  of  this  paper  under  the  side  condition,  that  the  antisymmetric 
part  of  the  solution  equals  Z(  in  equation  (14).  Surprisingly  enough  the 
answer  to  this  specialised  problem  is  easy  to  get.  By  elementary  operations 
one  finds  immediately  the  following 

Theorem  4  t  Assume  that  the  antisymsMtrlcal  part  of  the  circulation 

distribution  has  the  shape  of  a  lemniscate,  i.e.  that  Za  is  given  by 
equation  (14).  Then  there  is  exactly  one  circulation  distribution  with 
prescribed  lift  I,  and  roll  moment  R  for  which  the  induced  yaw  moment 
Y  vanishes  (no  adverse  yaw)  and  which  has  minimal  induced  drag  0.  It  is 
given  by  Zf  2  3  from  the  equation  (10).  The  corresponding  (minimal) 
Induced  drag  la  given  by  equation  (11)  . 
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4.  Bound!  For  The  Induced  Drag. 

In  this  Section  lower  and  upper  bounda  for  the  optima.  D  of  the  Induced 
ora*  o  are  given.  To  make  theae  r.  -alts  as  useful  as  possible  to  airplane 
engineer,  the  bounds  are  presented  for  the  optimal  drag  coefficient  cQ  as 
inequalities  of  the  form 


\  4  *  X  CD  ‘  "k  t0r  * 


2,  ...  • 


There  are  3  lower  and  4  upper  bounds  1,« 
defined  as  follows 


1-,  13  and  u,.  u2, 


u3,  u4*  They  are 


I  .=  c  2  +  8c  2  . 

I I  L  R 


l2  !’  "L 


+  3c  2  +  I  "Xcv  3cl°R  1  ' 


‘3  CL 


+  4cr2  ♦  HI  CR4  ♦  ||  (1 Ocy  - 


3clV 


2  , 


tr  Ac 


u  c  2  ♦  PC2  t  -  3cl)  for  cR  k  0. 

1  b  K  R 


t  SCR2  -  2  I  nicY  -  3  cL  cR  I  . 


0.2304  c 


«•  8  c  +  1-08 
L  R 


(1  + 


0. 16P  ,  .2 


(0.3  -  c, 


,2  >  <0'3  -  cr’ 


for  cy  -  0  and  cR  f  0.3  , 


2  2 

*  *  4  c 

r.  R 


♦-  /  1 6  c, 


R4  ♦  (»  X  cy  -  3  cLcR> 


22 


The  flret  bounds  11#  lj.  u^.  Uj,  Uj  are  sketched  in  Figure  9  for  Y  «  0 
together  with  the  solution  which  is  represented  by  a  dashed  line.  For  T  «  0 
the  whole  problem  is  homogenuous.  Therefore  «■  1  was  chosen  w.l.o.g.. 


Figure  9.  Upper  and  lower  bounds  u,,  u2.  Uj  and  1^,  lj  to  the  coefficient 
co  of  the  optimal  Induced  drag  as  functions  of  the  coefficient  cR  of  the 
roll  moment.  The  case  cL  s»  1,  cy  t-  0  is  treated.  The  dashed  line 
represents  . 


The  beet  bounds  which  were  found  in  this  paper  are  lj  and  u^.  They  have 
both  the  same  structure  and  differ  only  in  one  constant  (  24/25  in  lj  and  1  In 
u4  )  .  the  relative  error  between  them  is  less  than  2.1%  (l).  Hence,  for 
most  purposes  in  aircraft  engineering  these  easy-to-handle  bounds  are  com¬ 
pletely  satisfactory .  In  most  cases  it  will,  therefore,  not  be  necessary  to 
solve  the  problem  analytically  The  bounds  1^  and  u^  are  not  included  in 

Figure  9  since  they  are  both  too  close  to  the  solution  and  would,  therefore, 
interfere  with  the  picture. 


Proof  Of  The  Bound  Prop 


In  this  section  the  existence  theory  of  Section  5  will  not  be  used 

directly.  It  is  assumed,  however,  that  a  (not  necessarily  unique)  solution 

a  to  the  problem  of  this  paper  exists.  This  existence  is  nearly  trivial, 

since  D'(a)  is  a  continuous  functional  of  a,  bounded  by  zero  from  below.  For 

simplicity  the  notation  D  :«  D’(a)  with  the  corresponding  optimal  drag 

coefficient  c  is  introduced  . 

D 

Lower  Bounds 


This  bound  is  trivial  since 


I»au2  >  >  2?2  -  1,/sV. 

V*1 


Note  that  the  first  term  in  11  is  the  well  known  Prandtl  formula  (6)  for  the 


Ludwig  Prandtl  once  said:  "Es  gibt  nichts  Praktischeres  als  eine  gute 
Theorie"  l  (Translated:  "There  is  nothing  more  practical  than  a  good  theory") 


Ad  11  : 

D  :=• 


24 


induced  drag  of  an  elliptical  circulation  diatribution  (aee  Theoree  1).  It 


haa  been  uaed  by  airplane  engineera  for  aora  than  6S  yeara.  The  aecond  ten* 
in  11  ahowa  that  it  raquirea  at  leaat  an  induced  drag  coefficient  of 
8  cR2  to  produce  a  roll  aooant  coefficient  with  the  value  of  cR,  regardleaa 
of  the  value  of  cY-  The  bound  1^  correaponda  to  the  fonula  (8)  f roa  Theorem 
3.  The  value  of  t  X  equal*  1,  in  the  apeclal  caae  where  a^  -  L' , 
a2  -  R'.  In  thia  caae  Y’  -  3  L'R' . 

Ad  12  i  By  rewriting  O'  one  geta 

m 

4D'  i-  4  T  va  2 
l  v 

v-1 

m  m 

-  I  (  2v+1)  a  2  ♦  l  (2v-1)  a  2 

v-1  v-1 

-  a,2  ♦  I  (2v-M)(.v2  ♦  av+12)  . 

v-1 

By  adding 

t  2Y'  -  t  l  (2V+1 )  2  au  av+1 
v-1 

one  geta 

40’  ±  2Y'  - 


A  backtr ana formation  then  givaa  12 ■  Cbvioualy,  aa  can  be  aeen  from  Figure  9, 


•t  +  l  (2v+l ) (ay  t  .  )J 
v-l 


i,2  +  3  (a^  a2)2  . 
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For  large  value*  of 


the  bound  12  i*  better  than  1^  for  low  value*  of  cR. 
c p  the  opposite  la  true. 


Ad  1^  :  If  one  applies  the  Cauchy-Schwart  Inequality  to  the  Identity 


Y*  -  3  a,  a,  -  l  (2v+1)  a  a  , 
12  U  V  V+1 

v-2 


and  uses  the  inequality 


mm  <  a  ,v  ♦  „ 

v  6 


for  v  >  2 


one  get 3 


(V  -  3  a, a,)2  <  “  !»,  u2  I  (v+1)  a^ 


1  2 


v-2  v-2 

25  2  2  2 

—  (D* (a )  -  a1  )  (D* (a)  -  a1  -  2a2  ) 


By  solving  this  quadratic  inequality  for  D' (a)  -  a^  and  using  the  lower  bound 


1 ^  and  by  transforming  it  to  the  coefficients  one  gets  the  bound  1^. 


Tipper  Bounds 


To  get  upper  bounds  for  D  is  quite  easy.  Since  the  solution  a  is 
defined  by  a  minimum  problem,  any  approximation  a  to  a  which  satisfies  the 
side  conditions  (1),  (2)  and  (3)  gives  an  upper  bound  D'(a)  >  D  .  In  what 

follows  for  simplicity  the  sign  a  will  be  replaced  by  a* 


Ad  Uj  :  Let  a  be  defined  by 
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The  condition  (3)  than  reads 


V  -  3  a,  .2  ♦  5  «2  s3 
Hanes  for  R  ♦  0  one  gats 


*3  5  R'  "  5  L'  ‘ 


Inserting  thasa  coefficients  in  D'(a)  than  gives  u 


Ad  u2  t  Here  the  coefficients  of  a  are  dafinad  by 
av  0  for  v  +  J,  2,  V,  P+1  with  u  >  3  . 

The  condition  (3)  reduces  to 

V  -  3  .,.3  ♦  (211*1)  Vu+1  . 

hence  one  of  the  coefficients  a^  or  a^+^  rasalns  free.  By  chooaing  quite 
arbitarily 

^  \  +  ^  \+1  ■  0 
one  gets  then 


a 


M 


/  ATTT  1 3a  1  *2  rFT 

U  2u  •  1 


Inserting  these  coefficients  into  D'(a)  then  gives 


D' 


(a)  -  a,2  ♦  2  -j'  +  *  2  13  *1*2  '  Y>  1  ' 


where,  because  of  p  >  4  ,  the  coefficient  of  the  third  tare  is  bounded  by 


0.9938...  -  <  2/^T  <  1 

9  2  v  +  1 
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Raplacing  it  by  tha  upper  bound  1  and  transforming  back  than  ylelda  Uj. 


Ad  Uj  i  In  this  case  ay 


0  for  v  >  5  .  The  condition  (3)  then  gives 


a 


3 


Y'  -  3  a}a 
5  a^  +  7  a 


for  arbltary  a^  (provided  that  the  denominator  does  not  vanish).  The 
following  suggestion  for  a^  for  the  case  Y  “  0  has  been  given  to  me  by 
Mrs.  Norbert  in  Freiburg  i.Br./  GERMANY! 


* 


X  a . 


0.24  cr 

(0*3  -  v2 


for  cR  +  0.3  . 


Inserting  everything  in  D'(a)  then  gives  u^. 


Ad  u  :  For  finding  the  lower  bound  1,  the  Cauchy-Schwarr  inequality  was 

4  J 

used.  This  inequality  turns  into  an  identity  if  the  different  coefficients 
are  proportional  to  each  other.  This  leads  to  the  following  “Ansatr"  : 
Define  recursively 


av+1 


2v 

2v+ 1  av 


for  v  >  2  , 


where  the  constant  c  is  to  be  determined  later.  This  gives 


with 


cV”2  a2  for  v  >  2 

(2V-2) (2V-4)  ...  B  6  4 
(2V-1M2V-3)  ...  9  7  5  ' 


It  is  easily  proved  that 


With  the  aid  of  Stirling's  formula  it 
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for 


can  also  ba  shown  that  thara  exists  a  positive  constant  S  *  v  tj* 
all  v  >  2  .  Using  the  recursion  formula  once  more  gives 
<2v*1)  avsv+1  -  2  e  v  ay2  for  v  >  2 

and  therefore 


’  ’  3  *1*2  "  j  <2v+1)  *v*v*1  "  2  c  2  w  * 


v-2 


v-2 


-  2  c  D*  ( a)  -  2  a  c  l  v  6  ’  c 

2  v-2 


2  2v— 4 


for  any  real  value  of  V'  -  3a ^a2  there  is  axactly  one  solution  c  to  this 
equation  with  c2  <  1  .  This  cosms  from  the  fact  that 


r  ,  -  2  2V-4 

c  2.  v  o 

v-2 


is  a  monotone  function  of  c  which  goes  to  ±  •  aa  c  -♦  11.  This  can  be 

seen  from  the  inequalities 


e 


1  -  c  v-2 


1  -  c 


By  solving  the  quadratic  inequality 


If  -  3  .,.,1  <  4  a  2 

1  -  c 


one  gets  therefore  a  lower  bound  c*  for  |c|  with 


_1_  *_*/  *  /  4  a24  t  |y  -  3  a^l2 

Icl  c 


I  “  3  «,»2  I 
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By  inserting  this  in  the  shove  equation 


D‘  ( a) 


V  -  3  a1.2 


2  c 


and  transforming  it  back,  one  gets  finally  ■ 


Caeiparlson  between  and  1^ 


Using  the  abbreviations 


16  e 


(»  1  c„  -  3 


VV 


one  finds  for  the  relative  error 


U4  2  *3  ,  ✓  r  +  8  -  /  r  ♦  24  s/25 

*3  c  ^  ♦  4c  ^  +  S  r  +  24s/25 

L  R 

«  _ 8  (  1  -  24/25  )  _ 

/  r  +  24  s/25  (  /  r  ♦  s  +  f 


/  24/25  (  1  +  /  24/25  ) 


— - -  *  0.02062... 

2/6  (5  +  2/6  ) 


as  predicted  . 
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r  *  24  s/25  ) 


<  0.021  -  2.1  »  ] 

j 
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5.  Existence,  Unlquensas  And  Computation  Of  The  Solution. 

There  ere  several  approaches  to  a  proof  of  the  existence  and  uniqueness 
of  the  solution  to  the  problem  given  in  this  paper.  If  these  proofs  are 
constructive,  they  si so  permit  the  computation  of  a  solution.  In  the  paper 
[2],  the  author  used  the  theory  of  singular  integral  equations  to  prove 
existence  and  uniqueness  of  a  solution.  It  mas  even  possible  to  give  an 
explicit  expression  for  the  optimal  circulation  distribution  by  using 
Integrals.  Some  confutation  with  these  formulae  has  been  performed  in  a 
Master's  Thesis  ("Diplom-Arbeit")  at  the  University  of  Freiburg  l.Br. /Germany 
by  Mr.  Ramdane  Xedache  (unpublished).  Unfortunately  however,  the  numerical 
computation  of  these  integrals  (some  of  them  singular*  others  elliptic 
integrale  of  first,  second  and  third  kind)  was  very  etrenuous  work,  even  with 
today' e  feat  computers. 

Therefore,  an  algebraic  approach  is  used  in  the  present  paper.  The 
simplest  method  Is  to  replace  the  Infinite  series  In  the  equations  (3)  and  (4) 
by  finite  sums.  By  using  sufficiently  many  terms,  any  accuracy  desired  can  be 
achieved.  The  computations  for  the  results  presented  In  Section  3.4  have  been 
carried  out  by  this  method. 

In  what  follows,  however,  the  problem  Is  treated  in  a  different  way  by 
uelnq  the  "Lagrange  method*  2  It  will  be  shown  that  the  coefficients  a^  in 
a  solution  satisfy  a  certain  recurrence  relation,  with  this  result  their 
theoretical  structure  is  completely  understood.  What  remains,  is  to  determine 
a  certain  free  constant  in  that  relation.  This  can  be  computed  from  a  one- 
’ imenslonal  transcendental  equation. 
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5. 1  The  Lagrange  Method. 

In  the  problea  at  thla  paper  three  condition*  (1j,  (2)  and  (3)  have  to  be 
satisfied  with  arbitrary  constants  L,  R,  Y  or  L ' ,  R’,  V*  respectively.  The 
first  two  of  thee  are  linear  and  of  very  simple  structure.  To  solve  thee  gives 
no  difficulties  at  all.  Hence,  one  has  only  to  concentrate  on  the  one 
remaining  quadratic  condition  (3)  for  V'. 

In  the  Lagrange  method  the  following  fact  is  used:  Suppose  a  is  a 
solution  to  the  problem  (1)  to  (4)  .  Let  6  be  an  arbitrary  real  constant. 
Then  not  only 

r  2 

D'ta)  i-  l  v  a  •  minimum  for  a  ”  a, 
v-1  V 

but  also  the  same  is  true  for  the  following  quadratic  expression,  namely  that 

»  OO  ^ 

Q(a)  s«  j  v  ay2  +  8  J  (2v+1)  ayay+1  •  minimum  for  a  »  a  . 

v-1  v-1 

The  constant  8  is  called  the  Imgrange  parameter.  (Normally,  the  letter  X  is 
used/  but  in  the  present  paper  the  symbol  X  is  already  the  aspect  ratio.) 
Now,  by  rearranging  the  terms  in  Q  a  necessary  and  sufficient  condition  for 
a  is  reached  which  guarantees  that  a  is  a  solution. 

Let  a  e  A  be  a  solution  and  let  a  ■  a  +  b  £  A  be  an  arbitrary 
sequence  which  satisfies  the  conditions  (1),  (2),  and  (3).  Hence,  for  the 
first  two  coefficients  of  b  one  sees  immediately  that  b^  »  b^  ■  0.  By 
inserting  a  into  2  and  after  some  transformations  one  gets  the  following 
identity 
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0<a>  -  Q(a) 


P 


■  T.  ♦  T_  ♦  T. 


with  the  thra«  terms 


T,  1  \  l  B(2v+1 )  «v+i  ♦  2v  «v  ♦  B  (2v-1)  a^]  , 

v*3 

4  T  £  (2V  ♦  1)  (b  ♦  2  $  b  )2  . 

v-2 

m 

4  T3  i-  (1-4  a2)  I  <2v  ♦  1  )  b  2  . 

v-2  v ' 

Now,  a  la  a  solution  to  tha  problem  if  and  only  if  Q(a)  -  Q(a)  >  0 

for  any  a  c  A  satisfying  tha  conditions  (1),  (2)  and  (3)  and  that  Beans  for 

(almost)  any  b.  The  first  term  is  linear  in  b  .  It,  therefore,  has  to 

vanish  for  arbitrary  b  .  This  is  true  iff  all  the  square  brackets  are  aero. 
The  next  term  T?  is  quadratic  and  obviously  T2  >_  0  for  all  values  of  b  . 

The  same  is  true  for  the  third  term  Tj  iff  02  4  1/4.  If  one  puts  all  these 

facts  together  one  gets  finally  the 

Theorem  5  :  Necessary  and  sufficient  for  a  e  A  to  be  a  solution  of 

the  problem  (1),  (2),  (3)  and  (4)  are  the  following  three  conditional 


I)  a  satisfies  (1),  (2)  and  (3), 

II)  there  is  a  real  constant  B  with  I  B  I  <  1/2  such  that 
Hi)  a  satisfies  the  recurrence  relation 

(15)  B  (2v*1 )  ay+1  ♦  2v  av  +  8  (2v-1)  ay-1  -  0  for  v  -  3,  4,  ...  . 
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5. 2  Solving  The  Recurrence  Relation 

If  Be  (-1/2, 1/2]  and  if  aj  *•  R<  and  a^  is  an  arbitrary  raal 
constant,  then  (15)  can  be  solved  uniquely  by  recursion.  Hence,  there  is  no 
doubt  about  the  existence  of  a  solution  of  (15).  The  problem  in  question  in 
Theorem  5,  however,  is  different:  Is  it  true  that  for  any  values  of  L',  R* 
and  Y'  there  exists  a  constant  B  e  (-1/2, 1/2]  and  a  starting  value  a^  such 
that  the  sequence  a  generated  by  (15)  satisfies  the  conditions  (1),  (2)  and 
(3)  and  that,  moreover,  a  c  A  ?  If  yes,  are  the  constants  B  and  a^  then 
uniquely  determined  ?  It  will  be  shown  that  all  this  is  true.  But  it  is 
useful  first  to  treat 

The  Special  Case  B  »  0: 

In  this  case  the  recurrence  relation  (15)  has  for  all  starting  values 
a2  the  uniquely  determined  solution  *  0  for  v  >  3  .  Therefore  this  case 

is  equivalent  to  the  problem  treated  in  Section  3.2  with  the  solution  given  by 
Theorem  2.  Hence,  w.l.o.g.  B  f  0  can  be  assumed  in  what  follows. 

Obviously,  the  recurrence  relation  (15)  is  asymptotically  equal  to  the 
similar  relation 

(16)  B  ( 2  v  +  2 )  a  ♦  2v  a  ♦  S  (2v-2)  a  "0  for  v  3  3,  4,  ...  . 

v+1  v  v-1 

As  for  any  nondegenerate  second  order  linear  recursion  relation  both  (15)  and 
(16)  have  *•*■">  1  inear ily  independent  solutions.  The  two  solutions  of  (16)  for 
0  <  <  1/4  are  obviously 
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*u  aV/v  for  V  -  2,  3,  ...  , 

where  a  is  one  of  the  two  solutione  of  the  characteristic  equation 
2 

(17)  6  a  *  a  +  6  »  0  . 

These  solutions  are 

a,  (-1  +  /  1-402>/20  , 

and 

a2  (-1  -  / 1-4$2>/26  . 

Proa  these  expressions  ons  sees  easily  that  |a^|  <  1,  | |  >  1  for  f)2  <  1/4  . 
Therefore  only  the  solution  with  a ^  is  a  convergent  sequence,  moreover,  only 
this  one  solution  belongs  (  obviously)  to  the  class  of  admissible  sequences  K. 
Because  of  the  homogeneity  of  (16)  the  solution  wanted  is  then 

v— 1 

“l 

(18)  av  2  a2  — —  for  v  -  2,  3,  ...  . 

Inserting  this  into  the  condition  (3)  together  with  a^  “  X.'  and  a^  *  R' 
gives  the  equation 


Y' 


3  L'R' 


+  4  R 


.2 


m 

l 

v«2 


2v+1  2v-2 

v(v+1)  °i 


It  can  be  shown  rather  easily  that  there  is  a  unique  solution  a to  this 
equation  for  any  real  values  of  L',  R*  f  0  and  Y*.  Reinserting  this  solution 
into  the  characteristic  equation  (17)  then  gives  the  desired  value 
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6  !> 


1  +  a , 


If  R'  •  0  and  Y*  »  0,  then  obviously  the  Prandtl  solution  (5)  solves 
the  problem  of  this  paper.  Hence,  only  the  case  R'  «  0  and  Y '  f  0  remains 
to  be  treated.  This  will  not  be  done  here. 

Up  to  this  point  only  the  asymptotic  recurrence  relation  (16)  has  been 
treated.  It  can  be  solved  explicitly  by  (18).  It  can  be  shown  that  (18)  is 
also  an  asymptotic  solution  to  the  original  equation  (15).  This  will  be  done 
elsewhere.  Hence  one  gets  as  a  result  the 


Theorem  6;  For  any  combination  of  real  numbers  L‘,  R‘  t  0  and  Y'  there  is 
exactly  one  solution  a  €  k,  8  to  the  recurrence  equation  (15)  such 
that  the  conditions  i)  and  ii)  of  Theorem  5  are  satisfied. 


5.3  Uniqueness 

It  can  be  shown  that  there  is  a  solution  to  the  problem  of  this  paper 
also  in  the  case  R  «  0,  Y  f  0,  which  was  not  treated  above.  But  in  this  case 
the  solution  Is  not  unique;  in  opposition  to  the  above  result.  This  (.an  bee 
seen  very  simply  as  follows: 

-  •  a  “  {  a.  ,  0,  a. ,  a.  ,  . . .  }  be  a  solution  to  the  problem  (1), 

1  3  4 

(2),  (3)  and  (4)  with  R'  -  0,  l.e.  a2  =  0  .  Define 
a*  :»  {  a1#  0,  -a3  ,  ...  }  . 

Then  a*  obviously  also  satisfies  the  conditions  (1),  (2)  and  (3)  .  Moreover 
D' (a)  “  D'(a*)  which  shows  that  a*  is  a  sc  "ton  of  the  problem,  too.  It  can 
be  shown,  moreover,  that  these  are  the  only  -wc  solutions,  i.e.  that  no  other 
solutions  exist.  The  proof  of  this  will  not  b*-  gi  'er>  here. 
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5.4  Generating  Function 


(Thi*  part  has  been  proposed  to  me  by  Professor  Ben  Noble,  former 
director  of  the  Mathematics  Research  Center.) 

Let  z  be  a  complex  variable  and  define  the  complex  function  f 
power  series 


f  ( z ) 


v-2 


2v 


Such  an  fuction  f  is  usually  called  a  ’generating  function"  of  the 
sequence  a  .  Then  inside  the  disc  of  convergence  one  gets  : 

OD 

(f(z)/z)'  -  l  ( 2v+1 )  av+1  z2J  , 
v-1 

oo 

*  f  (z)  »  l  2  v  a  *2V  . 

v«2 

o° 

Z2(Z  f ( z) ) •  -  l  (2v-t)  av_1  z2U  . 
v-3 

Assume  that  the  coefficients  av  satisfy  the  recurrence  relation  (15) 
by  the  above  three  equations  it  follows  immediately  that  the  function 
satisfies  the  following  linear  differential  equation 

z(B  +  Z2  ♦  8z4)f'(z)  +  B  ( z  4- 1 )  f  ( z )  -  (az(3B  +  4z2)  +  5Ba3z2)z4  . 


With  the  Initial  condition  f ( 0 )  •  0  it  has  the  unique  solution 


f(z) 


z 


-H(t) 

e 


k(t)  dt 


by  the 


.  Then 
f 
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with 

h<z)  : * 

H ( z )  := 

Mz!  :  = 

From  this  generating  function  the  results  of  Theorem  6  can  also  be  derived. 
This  will  not  be  shown  here. 
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a2z(3B  +  4z  )+  5  3  a^z 

2  4 

8  ♦  z  +  6  z 
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